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Abstract 

Position distributions of constituent particles of the perfect Bose-gas trapped 
in exponentially and polynomially anisotropic boxes are investigated by means 
of the boson random point fields (processes) and by the spatial random distribu- 
tion of particle density. Our results include the case of generalised Bose-Einstein 
Condensation. 

For exponentially anisotropic quasi two-dimensional system (SLAB), we obtain 
three qualitatively different particle density distributions. They correspond to the 
normal phase, the quasi- condensate phase (type III generalised condensation) and 
to the phase when the type III and the type I Bose condensations coexist. 

An interesting feature is manifested by the type II generalised condensation in 
one-directional polynomially anisotropic system (BEAM). In this case the particle 
density distribution rests truly random even in the macroscopic scaling limit. 
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1 Introduction and Main Results 

The microscopic position distribution of constituent particles of quantum gases in sta- 
tistical mechanics can be described by Random Point Fields (RPFs, or Random Point 
Processes). For bosons this idea goes back at least to [F l IFF) . Determinantal and per- 
manental processes describe, respectively, the position distributions of ideal Fermi-and 
Bose-gases [M751 Ly . It is known that position distribution of bosons in the state of 



Bose-Einstein Condensation (BEC) can be described by the corresponding permanental 
RPF explicitly |TIbt lEf ITZaj . For a introduction into the theory of general RPFs as 
well as the fermion- and the boson- RPF and their applications in Statistical Mechanics, 
see, e.g., [DVl iGioYl ISTl HHUfc]. 

It is also possible to prove the Large Deviation Principle (LDP) for the boson RPFs 
in the condensed and non-condensed phases \ST\ ITZb] , extending some known results in 
this field [BruZj . although the critical point needs a separate analysis. The proof of the 
LDP involves large-scale test functions for the generating functional of the boson RPF, 
which one must scale properly for different phases, see discussion in [TZb] . Alternatively, 
the large-scale behavior of position distribution, or macroscopic (or mesoscopic) position 
distributions, can be treated as a particle density, which can be described by Random 
Fields (RFs). For example, RF which represents macroscopic boson density is obtained 
by scaling the RPF for boson gas trapped by external potential, see |TZa[ IT]. However, 
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in the last case the RF yields a deterministic density distribution in the semiclassical 
limit, which is similar to the law of large numbers |TZb] . 

Recently, the study of boson gases trapped in quasi one-dimensional and quasi two- 
dimensional regions have been developed |BZ] with the aim to elucidate physical argu- 
ments and experiments, which indicate the existence of generalised (gBEC), or quasi- 
BEC (qBEC), boson gases trapped in extremely anisotropic systems |KvDl IvDKl IGAl 
IHD] . The unusual properties of highly anisotropic boson systems are known since eight- 
ies and they can be related to the notion of generalized BEC by van den Berg-Lewis- 
Pule (vdBLP-gBEC) [BLP] . Although their observat ion and classification of the gBEC 
in types I, II and III seems purely mathematical at the first glance, it is accepted to be 
quite relevant to the interpretation of the experimental data in extremely anisotropic 
systems, such as SLABs, BEAMs and CIGARs, see discussion in |BZj and the recent 
paper |MuSaj . 

However the case of the anisotropic CIGAR shape, which is the most interesting for 
traps, is not covered by the van den Berg-Lewis-Pule theory. For the first time this 
case was rigorously studied in |BZj . It gives a qualitative description of the experiments 
|vDKt \QA\ IHD] . including the existence of two critical points and transition between 
type-I and type-Ill condensates, i.e. between normal conventional BEC and "quasi- 
condensation" . The peculiar property of the CIGAR systems seems to come from their 
exponential one-direction anisotropy. Although the BEAM anisotropy may produce the 
type III condensation, it does not split the single critical point into two points, if one 
does not modify external potential in the anisotropic direction |BZ] . 

The aim of the present paper is to study the features of the generalised BEC phases 
from the view point of position distributions. To this end we consider RFs which describe 
scaled particle density as well as RPF for BEC in two cases: the exponentially anisotropic 
SLAB and algebraically anisotropic BEAM which are obtained by the thermodynamic 
limit of Casimir prisms. 

As mentioned above, the SLAB system manifests three phases: normal phase (with- 
out BEC), condensation of type-Ill and coexistence of condensations type-I and type- 
Ill, separated by two critical points. So, first we construct the RPF, which describes 
the microscopic position distribution of constituent bosons, for these three domains. 
This result clearly demonstrate the difference between the normal and the BEC phases. 
However, the qualitative difference between two BEC phases (type-Ill and the mixture 
of type-Ill and type-I) becomes evident only after contraction to the particle density, 
obtained RPs by scaling the boson RPF. We find that the density distributions are 
essentially deterministic, i.e., the corresponding RFs are non-random. This result is 
similar to the behaviour of the particle density distribution in the "semi-classical" limit 
for the boson mean- field model in a weak harmonic trap jT]. 

In contrast to the SLAB model, the BEAM is an algebraically anisotropic case of 
Casimir prisms. A specific choice of anisotropy rate gives BEAM a stronger randomness 
to the particle density distribution than the SLAB. It is known that anisotropic BEAM 
has only one critical point |BLP| IBZj . Our specific choice of the anisotropy rate (see 
Section 1.3) implies that this point separates the normal phase and the type-II conden- 
sation. Since the type-II condensation is in sense a type-I condensation smeared out 
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over a countable infinite number of states |BLPl IZBj . it is expected that the particle 
number fluctuations in type-II condensate should be higher than the type-Ill conden- 
sate phase. jBLPj IBLL] Our result shows that for this kind of BEAM, the macroscopic 
particle density RF obtained by scaling (similar to the SLAB case) from RPF gives 
non-trivial fluctuations described by squared gaussian distribution. 

In the next subsection, a brief introduction to boson RPF is presented in order to 
describe the position distribution of the free boson gas in three dimensional anisotropic 
boxes. In subsection 1.2, we formulate theorems on RPFs and RFs concerning micro-, 
macro- and mesoscopic behaviour of position distribution of the boson gas in exponen- 
tially anisotropic SLAB as well as their corollaries. In subsection 1.3, we collect the 
corresponding results about algebraically anisotropic (in one direction) Casimir prisms, 
that for a spacial rate converges to the BEAM. Section 2 and 3 are devoted to proofs of 
the theorems for SLAB and for BEAM, respectively. 

1.1 Preliminaries: Notations and Definitions 

(a) Random Point Fields (Processes). Let £ be a locally compact metric space 
serving as the state-space of the point configurations £ C E. By 03 we denote the 
corresponding Borel a-algebra on E and by 03o Q 53 the relatively compact Borel sets 
in E. We denote by /i a diffusive (i.e. /i(x) = for any one-element subset x G E) 
locally finite reference measure on (E, 03). (The standard example is the Lebesgue 
measure fi(dx) = dx on (E = M d , 03).) 

We denote by Qe the subspace of locally-finite point configurations {£ C E} : 

Q E ■= {£ c E : card(£ fl A) < oo for all A G 03 } . 

Hence, for any A G 03 one can define a subspace of the point configurations Qa : = {£ £ 
Qe :(C A} and the mapping n\ : £ \-> £ fl A for the corresponding projection from Q E 
onto Qa- Then counting function: N\ : £ H- card(7TA(£)) is finite for any A G 03o- 

Now one can introduce the notion of the spatial random point field (RPF) on M. d as 
locally finite discrete random sets £ C M. d , i.e. such that N\(£) < oo for A G 03o. Since 
below we use the Laplace transformation for characterisation of the RPFs, we need a 
more elaborated general setting. 

Let 5 X denote the atomic measure on 03 supported at one-element subset x G E. 
Then any configuration of points £ G Qe can be identified with the non-negative integer- 
valued Radon measure: A^(-) := J2{xe^} ^(*) on ^ ne Borel a-algebra 03. Hence, \%(D) = 
Ne>(£) is the number of points that fall into the set D G 03o for the locally finite point 
configuration £ G Qe- Recall that Cq(E)*, which is dual to the space of continuous 
on E functions Cq(E) vanishing at infinity and equipped with the uniform norm, is 
isometric (by the Riesz representation theorem) to the space Ai(E) of Radon measures 
on E. By this isometry the weak-* topology on Cq(E)* yields the vague topology on 
M.(E). Then identification of M.(E) with the set of Radon measures induces on the 
point configuration space a topology, turning Qe into a locally compact separable 
metric space with the corresponding Borel a- algebra 03 (Qe)- Note that if $ is the 



4 



smallest a-algebra on such that the mappings are measurable for all A G 23o! 
then £ = *8(Q E ), see e.g. [DV] . 

Definition 1.1 A random point field (process) is a triplet (Qe^{Qe)^ v ), where v is 
a probability measure on (Qe,%$(Qe)) ■ Us marginal on Qa is defined by the probability 
measure := v o n^ 1 . 

Note that the process defined above is simple, i.e. the random measure almost 
surely assigns measure < 1 to singletons. 

(b) Correlation Functions and Laplace Transformation. For the marginal mea- 
sure is\ we consider the Janossy probability densities {ja jS (xi, . . . ,x s )} s > , |DV] . Here 
j A )S= o(0) = ^a({£ : -^a(0 = 0}) and for s > 1 it is a joint probability distribution that 
there are exactly s points in A, each located in the vicinity of the one of x\, . . . , x s , and 
no points elsewhere. By construction the Janossy probability densities are symmetric 
and verify the normalization condition 

°° 1 f 

5^— / /i(dxi) . . .fJ,(dx s ) jA jS (x!, . . . ,x s ) = 1 , (1.1) 

with a standard convention for s = 0. Then for any measurable function F on Q\ with 
components {F s } s > one gets |DVj : 

/ VA(d£)F(£) = y~] — I ii(dx 1 ).../i(dx s ) j AjS (xi,...,x s ) F s (x!,...,x s ) . (1.2) 
JQa s=0 s[ Ja s 

These joint probability distributions (correlation functions) serve for a very useful 
characterization of RPFs by Laplace transformation. Let / : E — >■ K + , be non- negative 
continuous function with compact support. For each / one can define by 

(/,£):= / A*(<fa)/(s) = , (1.3) 

the measurable function: £ i-» (/, £) on Q^. Then by virtue of (11. 2p the Laplace 
transformation of the measure v\ for a given / takes the form 



E„ A (e-<«>) = / i/ A (de) = (1-4) 

^ — / /i(cixi) . . .fi(dx s ) j\, s (xi, . . . ,x s ) ~[e~ f{Xj) . 
s =o s ' ^ As j=l 

The most fundamental example of RPF is the Poisson point process ir z (d£) on E = W 1 
with Lebesgue measure fi(dx) = dx and the intensity function z(x) > 0. For this RPF 
its marginal on Qa is defined by the Janossy probability densities: 

s 

{j A , a (xi,...,x s ) =e~JA Y[ z { Xj )} a >o . 

3=1 
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Then taking into account fll.4p one gets for any non-negative continuous function / 
with compact support the corresponding Laplace transformation (generating functional) 
expressed by the well-known formula: 

En^e-UM) = exp l- J dxz(x)(l - e~ fix) )\ , (1.5) 
for extension to infinite configurations Q-^d. 

(c) Permanental (boson) RPF. In the box A = ]T^ =1 [— Lj/2, Lj/2] C M 3 , we consider 
a quantum mechanical system of identical free bosons of the mass m with the one-particle 
Hamiltonian 

H A = -^A D , (1.6) 
2m 

where denotes the Laplacian operator with Dirichlet boundary conditions in the 
space L 2 (A) . The spectrum and the eigenfunctions of this self-adjoint operator are 
given by 

e = ^ /7rfcj \ 2 

k 9m 2-^< \ T. . 
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k = (ki, k 2 , k 3 ) G N 3 j , (1.7) 



keN > . (U 



Then the corresponding trace-class valued (Gibbs) semigroup Zag| and its kernel are 



(G A (/3) := e"^ A e Ci(L 3 (A))}^ , G A (frx,y) := ]T e-^0 k , A (x)0 k , A ( 2 /) . (1.9) 

keN 3 

The random point field defined by this kernel is the Boson (BRPF), i.e. the perma- 
nental RPF |DV] . It describes, in correspondence with the RPF paradigm, the position 
distribution of non-interacting bosons in the box A. If the Bose-gas is in the grand- 
canonical equilibrium state at fixed temperature and chemical potential fi, then the 
BRPF generating functional (cf. fll.5p ) has the form [TZaj . [T] : 

JQa -AflW 



Here denominator is the perfect Bose-gas grand- canonical partition function for \x < 0: 

™ r e n, 

n=0 ' 

Det [1 - e^G A ] _1 , (1.1F 



,(/i) := ) I dx 1 ...dx n — j-Per{G , A (/3;Xi,a; i )}i^ iKr 



and the numerator 

j(/i) ■= ^2 dx x ...dx n — j-e" E "= l/(Xs) Per{G A (/3;x i ,x i )}i s ;i iKn 

n=0 ^ A ™ U ' 



Det [1 - e^e^GA]- 1 = Det [1 - e^G A /2 e~ f G 1 / 2 }" 1 , (1-12) 
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where Det denotes the Fredholm determinant and Per stands for the permanent of the 
corresponding operators. The both representations (11. lip and f ll . 1 2j) via determinants 
are due to Vere- Jones' formula generalized to the trace-class operators G A in Theorem 
2.4, [ST] . 

Notice that in the grand- canonical ensemble the BRPF configurations Q A are infinite. 
To consider in A the finite point BRPF ^a,«(c?0 one has to consider the canonical 
ensemble |TIa| ITIb] . Then the probability density of the n-point configurations in A is 
defined by 

PA, n (xi, ■ ■ ■ ,x n ) = —— — ^Per{G A (/3;£j,£j)}i^ i: Kn , (1-13) 

where 

Z Afl (n) := Tr Hn e ~ pH ^ n = [ dx x ...dx n — Per{G A (/3; x u x i )} 1<i)Kn , (1.14) 

J A« n - 

is the n-boson canonical partition function defined by the trace over the Hilbert space 
1-L n := L 2 ym (A n ) of symmetrized function |TIaj . Here the sum 

h 2 

l<j<n 

defines the n-particle Hamiltonian in the space L 2 (A n ), cf. (II. 6ft . Similarly to (II. 12ft we 
also define 

z A,f(n) := / dx 1 ...dx n — e~ T < n ^ f{Xs) ¥eT{G h (f3]x u x j )} 1 ^^ n . (1.15) 

Then according to construction presented in (a), the n-point BRPF /XA,n(^0 i n A 
is induced by the map (xx, . . . ,x n ) >->■ ^™ =1 (•) £ Qa,u and probability measure on 
A n with density (11. 13ft . The BRPF generating functional in canonical ensemble has the 
form [TIairnia] : 

Z Kn [f] = [ wUdO e-™ = §^7 > (1-16) 

Note that by (II. 2p and (11.41) the Janossy probability densities of the finite-volume 
BRPF for a fixed grand-canonical temperature (3 and the chemical potential fi are 

jA )n (P,fi]x 1 ,...,x n ) := - -. e n ^ Per{G A (/3; x h Xj)} 1<iJ<n , (1.17) 

see (ll.lUp . Here {n = 0, 1, . . .} and Per{GA(/3; 2 ? j)}i<tj^n=o = 1- 

Since the function / is supposed to be non-negative continuous with a compact 
support, the fundamental properties of the (Fredholm) determinant and (ll.lOp imply 
that: 

Zl C Jf] = Det + e~fK A (zWl - e"/]" 1 , z := ef*> . (1.18) 
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Here for fi < 0, the positive operator K\{z) := zG\{l — zG A ) 1 > 0, is self-adjoint and 
belongs to the trace-class Ci(L 2 (A)) of operators on L 2 (A). We also denote by 

K AJ (z) := Vl-e-fK A (z)y/l-e-f E Ci(L 2 (A)) , (1.19) 

the trace-class operator in determinant (I1.18p . 

In order to make the generating functional Z 9 A Af\ well-defined, i.e. to ensure the 
convergence in the numerator and the denominator [stability), the chemical potential 
/i must satisfy the condition e^^Gx < 1, i.e.,/i < €\ — Y^j=i K 2 h 2 /2mLj 2 . Here 1 = 
(l,l,l)eN 3 . 

Hence, from now on we put fi := ewi^i) — A, where for stability we assume that 
A > 0. Since below we examine the behaviour of BEC in the thermodynamic limit 
{L — > oo : Li,L 2 ,L 3 — > oo} for anisotropic vessels, we assume that A = Ax, may 
depend on L±, L 2 , L% and that the value of (LiL 2 L 3 ) A^ is bounded. The latter ensures 
boundedness of the total particle density: 



P 



V ^ , (1.20) 



LlL ^ » ex V P[(hH 2 /2m) E; =1 (*? - 1)/LJ + A L 



in the grand-canonical ensemble (/3,/x), including the limit L — > oo. Here {A^ = 
Al(/5,p)}l is solution of the equation (I1.20I) . [ZB\, for a fixed bounded density p. 

Notice that the integral kernel of the operator K A in the grand-canonical ensemble 
(/3, fi = e\ — A) is given by 



K A (x, y) := Y M*)Mv) (1 21) 

k ^a exp P[(h^/2m) Ej =1 (k] - l)/L) + A] - 1 

Since the series in the right-hand side of f II . 2 1 j) is uniformly convergent, the kernel is 
continuous. 

(d) Equivalence of ensembles and the Kac probability distribution. To make a 
contact between canonical (/3,p) grand-canonical (/3,/i) ensembles we use to the gener- 
ating functional f ll.lOp for the case of f(x) := £Ia(;e)/|A|, t > 0. Here Ia{x) is indicator 
of domain A with volume which is V := |A|. Then by virtue of (jl.3p the functional 
(II.IOP takes the form 

-/3(e k -/i) 



(e -t» A ®/V) = ZA(P,H-t/(0V)) _ -,-r 1 - 



where E / 3 iAt (-) is expectation in the Gibbs grand-canonical ensemble (/3, /i) and Sa(/3, p) = 
"A,o(/ i ) (II- lip . Notice that fll.22p has also representation: 

^A*- tNm,v ) = e ^i^r e ~ WF ' (L23) 



n=0 



where Z\(0,n) = Za,o(^) is defined by (j!.14j) . 



S 



For domain out of condensation, i.e. for p < 0, the thermodynamic limit of (I1.23P 
follows directly from representation (11.221) : 

POO 

Km E^e-^COA') = / dx5 p(M (x) e~ tx = e~^ M . (1.24) 

V->oo J Q 

Here 

1 f d 3 k 

PiP^) := (27Fy R3 e^-l' (L25) 

is the particle grand-canonical density for p < 0, which reaches its maximal (critical) 
value at p = 0: 

p c (/3) :=p(/3,p = 0) (1.26) 

Note that (ll.24p expresses the strong equivalence of ensembles |ZPj for the perfect Bose- 
gas in this one phase regime: 

POO 

= / dx5 p (p^)(x) Ep iX (A) = Ep )P ^^ (A) , (1.27) 

where Ep tX (-) is the canonical quantum Gibbs state for the total particle density x, and 
A, A are corresponding observable in grand- and in canonical ensembles |BNZj . 

Hence, for p < p c (/3) the solution of equation (ll.20p gives Pl(/3, P) such that the limit 
liim/^oo/izXAp) = M/3,p) < °) and p(/3,KP,p)) = P b y (OS]). 

Definition 1.2 We define the chemical potential ~p := p(/3,p) as the inverse function 
of density \1.25\) , \1.2b]) , with extension: p(/3,p > p c (/3)) := 0. Hence, p(/3, p(/3, p)) = p 
and p(/3,p(/3,p > p c (P))) = pdp). 

The phase with condensation corresponds to p > p c (P), i-e. limy^oo p L (/3, p) = 0. 
Now to deduce the asymptotics of the solution p^(/3,p) one has to separate from the 
sum (11 . 201) the most singular term 1 = (1,1,1). Let Pl(/3, p) := t\ — A^(/3, p). Then for 
thermodynamic limit (in the Fisher sense [Ru] ) it has the form: 

Mp)^^- ^_l mv + oiy^). (1.28) 

The asymptotics fll.28[) ensures that the single term 1 = (1, 1, 1) gives particle density 
Po(P) '■= P — Pc(P) > (Bose-Einstein condensation), whereas the rest of the sum 
converges to the integral fll .25[) for p = 0, i.e. to p c (/3). Therefore, by (I1.23P and fll .28[) 
one obtains the limit 

Km E MM (e- tN ^y y ) = / dx —— e" te (1.29) 



where limy-^ Pl(/3, p) = 0. 
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Remark 1.3 Representations \1.21$ , U.29\) and Definition \1.2\ motivate the notion of 
the Kac probability kernels {JCp^tp^ix, p)} p >o, see ILePuf : 

fc/3,jc(/3, P )(x,p) ■= Q+{pc{P) ~ p) Sp^ix) + (1.30) 

9_{x-p c (J3)) e -(-p c (/3))/(P-P C W) 
0-(p-Pc(P)) 7^ ■ 

p- pM 

Here 6 + (z > 0) = 1, 6+(z < 0) = and 9-{z > 0) = 1, 9_(z < 0) = 0. 

Application of (jl.3Uj) to calculation of expectations in the grand- canonical ensemble 

poo 

^p,v=np,p)( A ) = dxfCp tlx ^ t p)(x,p) Rp >x {A) , (1.31) 
Jo 

shows that when ~p(/3,p > p c (/3)) = 0, i.e. in the phase with condensate density po(/3) = 
P ~ Pc(P), the grand-canonical Gibbs state 

POO 

®p,ji(p,p>p c (P))=o(A) = / dxKp >0 {x,p)Ep tX (A) , (1.32) 

Jo 

is convex combination of the canonical states {W,p jX (A)} x>p rpy Here the measure 
dx )Cp s o(x, p > p c (/3)) defines the Kac probability distribution with support (p c ((3), +oo), 
see fOCTjl . 

Remark 1.4 Since the grand-canonical particle density N\(^)/V is non-negative ran- 
dom variable, the Laplace transforms \1.24\ ) and U.29\) defines characteristic function 
of the Kac distribution IfFM) . TBTPj . Notice that e -«*p(^./*) j ^ < o, and {l + t(p- 
Pc(P))}~ a , P > Pc(P) are characteristic functions for any a > 0. Therefore, two proba- 
bility distribution with densities: 

fcp,JZ(P,p)( X 'P) = 6 p(/B,Jl)=p( X ) °+( X ) > P < PM > ( L33 ) 

see Definition \1.2\ and 

x /(p- Pc (l3)) 

/Cfjg^p) := 9_(x) , p > ptf) , (1.34) 

P Pc{P ) 

are infinitely divisible fSB/ . 



Notice that (11.331) coincides with the Kac probability kernels fll.30p for p < p c (/3). On 
the other hand, by virtue of f 1 1 . 3 j) and (jl.34p the Kac probability kernel for p > p c (/3) 
is convolution of two probability measures: 

fcp,iz(p,p>pcm( x > P) = (^MP) * ' C iW 1 Pc))(s, P) = (1-35) 

d v tfM^ ~ y > K p,w,pc)(yi p)= I d v zfLwfc ~ y > p ) 5 pc{p)(y) > 
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see Definition 11.21 Since the convolution of two infinitely divisible distribution is again 
infinitely divisible, we deduce this property for the Kac distribution (jl.30p also in the 
BEC phase with condensate density po((3) = p — p c ((3) > 0. 

(e) Spacial distribution of Boson particles. Since our aim is to analyse of the 
space distribution of boson for different initial geometry, we would like to indicate how 
sensitive might be the Bose condensate shape as a function of the boundary conditions 
in the initial restricted geometry. This implies certain specificity of position distribution 
of bosons even in the macroscopic scale and in the thermodynamic limit see e.g. |VVZ] . 
(Verj . 

Recall that for the periodic boundary conditions the spectrum and the eigenf unctions 
are different to (11.7j) and (11.81) for Dirichlet boundary conditions: 



s k = — V — -i k=(k 1 ,k 2 ,k 3 )ez i , (1.36) 




^,a(z) = J— exp [i-r^-ixj + Lj/2) k G Z > . (1.37) 



Since the torus is wrapped, the kernel of the corresponding operator K\ is translation- 
invariant: 

1 TT 3 e i27rk j ( x j-yj)/ L j 

K per (x - v) ■= j=1 (1 38) 

A ^ exp P[(h*(2K)*/2m) £j =1 k V L ) + A l " 1 ' 

where A = —p > 0, since in this case £k=(o,o,o) = 0. 

Remark 1.5 The kernels M.21\) and U.38\) produce different "global" position distribu- 
tions of bosons and the shape of the condensate. It is visible from ground state density 

8 3 

|0k=l,A(^)| 2 = 1-1-1- I [ COS 2 (7TX,/%) , (1.39) 

L1L2L13 

for Dirichlet boundary conditions ( fi.gj) and 

i^ k =i,A(x)r=-^— , (1.40) 

for the periodic boundary conditions. Hence one has to anticipate the difference in the 
local (in the vicinity of x = 0) as well as in the scaled space distributions. Below we 
essentially stick to the Dirichlet boundary conditions. 

(f) BRPF theory of the conventional Bose-Einstein Condensation | ITIa|, ITlb] . 

A peculiarity of the Tamura-Ito approach |TIa| ITlb] is that BRPF is studied in the 
canonical Gibbs ensemble f II . 1 6 j) . To take into account the phenomenon of the Bose- 
Einstein condensation (BEC) they consider thermodynamic limit of the BRPF, when 
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A — > M rf for a fixed particle density n/\A\ = p. The BRPF is completely characterised 
by the limit of generating functional (I1.16p . and in |TIa] . [Tib] this limit is established 
along the family of A — > ~R d , which are non-anisotropic. This is a standard precaution 
to reach the conventional bulk properties and in particular the conventional BEC. It is 
known as the Fisher limit, for sufficient conditions see e.g. [Ruj . 

The first result jTIa] concerns the case of subcritical density of bosons p < p c ((3), 
cf.(d). Then for any non- negative continuous function / with a compact support the 
finite volume canonical BRPFs {/iA,n.}A converge weakly to the random point field fi p 
such that the limit generating functional f 1 1 . 1 6 j) takes the form: 

Z p [f] : = lira Z A>n [f] = [ (*„(<%) e"<«> = Det [1 + K f {z{p))]- 1 . (1.41) 



Here K f (z) : = Vl - e^K{z)\fl - e~f and K(z) : = zG p (l ~ zGp)' 1 , where Gp := e^ A 
is the heat semigroup in the whole space M. d , cf. (11.191) . and z(p) is a unique solution of 
equation 

1 f d 3 k 

P ~j2^L z-ieP«<-l ' (L42) 

for the total particle density p < p c (/3). Similar to the grand-canonical ensemble, see 
fll.25p . fll.26p . \im p ^ Pc (i3) z(p) = 1 and equation (11.421) has no solution for p > p c (f3)- 
The excess of the particle density po((3) := p — p c {fi) > accumulates in the ground 
state 1 = (1, 1, 1) (Bose-Einstein condensation), although absence of the explicit formula 
(ll.20p makes this analysis more complicated than in grand-canonical approach [PuZ] . 

In spite of these difficulties Tamura and Ito succeeded to find the limit of the canon- 
ical BRPFs {pa,u}a in the phase with condensation |TIb] : 

Z p [f] := lim Z A , n [f\ = [ ^ = ( L43 ) 



exp{-(p - p c (/3))(yT^7, [1 + Kjjl^VT^f)} 
Det [1 + K f (l)} 

Their detailed analysis of the canonical limit shows that the structure of generating 
functional (I1.43P is due to convolution: pf p = p Pc (/3) * A t p_ Pc ( J g) of two boson random 
point processes. The first one is (11.411) for the critical density p = p c ((3), whereas the 
second one ^° po ^ corresponds to the BRPF of particles in the condensate with density 
p ((3) :=p-p c ((3)>0: 

f e-<^> = exp{-r( v / l-e-/. \ Ji- e -f)} , r > 0. (1.44) 

Finally, taking into account (ll.lOp and definitions (Il.lip . (ll.l2p together with (I1.14p . 
(I1.15p . (jl.l6p . we obtain representation 

2Uf] = £ 2aM . (L45) 

n=0 
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Notice that discrete measures in (11.231) and (I1.45P are identical. Therefore, in the ther- 
modynamic limit the grand-canonical generating functional Z g ^ = -^ s [f] is related to the 
canonical one by the Kac transformation: 

/■oo 

Z »=HMW= dxK f!>mp) (x,p)Z x [f). (1.46) 

J 

Then taking into account (jl.30p this implies 

Z°t mp) [f] = 8 + ( Pc (f3) - p) Det [1 + KMpW 1 + (1-47) 



l + (p- Pc (/3))(v / T^7, (1 + K f (l))^VT^f) ' 



Generating functional (I1.47P expresses the Tamura-Ito results |TIa[ ITIb] in the grand- 
canonical ensemble. 



1.2 Microscopic and macroscopic distributions of constituent 
bosons in exponentially anisotropic SLAB 

(a) Two critical points and generalised BEC in anisotropic SLAB. First, we 
describe the space distribution of the constituent bosons in the thermodynamic limit 
of exponentially anisotropic (in x\, ^-directions) SLAB with Dirichlet boundary condi- 
tions. 

To this end we set L\ = L2 = L e aL , L3 = L for a > and consider the limit L — > 00. 
We also denote this parallelepiped by A = and by Vl = L 3 e aL its volume. 
Let p < 0. Then for the limits of the Darboux-Riemann sums fll.2Qj) one gets: 



lim 



1 E 



Vl t£ exp(3[(hV/2m) - 1)/L] + A] - 1 



lim — 



1 



L ^ Vl exp PW^/^m) e; =1 (k] - + a] - 1 

H 3 k 

■=: P (P,p = -A) . (1.48) 



Here p(/3, p) is the the grand-canonical particle density mean- value with finite critical 
density p c ((3) : = p(ft,p = 0). Since (jl.58p yields: 

lim V \ = p c ((3) , (1.49) 

^°°V(M 2 ,i) exp/3[(fr 2 7r 2 /2m) Y? j= i {k] - 1)/L?] - 1 

by virtue of ( 11 ,20j) for p > p c (j3) we obtain that: 

T,— inn V t < • 



tl iwrL e^y^M) - 1 = &, Ajz ln[/3AL(/3 ' p)] " ■ (L50) 
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This implies the asymptotics: 

A L ((3, p) = i e - x >- p < mL + o(e~ AL ) . (1.51) 

Notice that representation of the limit (ll.50p by the integral is valid only when A|(p — 
p c {/3)) < 2a. For p larger than the second critical density |BZj : 

p m {0) ■= Pc{P) + 2a/Aj (1.52) 

the correction Al(P, p) must converge to zero faster than e~ 2aL . Now to keep the differ- 
ence positive: p — p m ((3) > 0, we have to return back to the original sum representation 
(ll.49p and (as for the standard BEC) to take into account the impact of the ground state 
occupation density together with a saturated non-ground state (i.e. generalised) con- 
densation p m (/3) — p c (j3) as in (11.501) . For this case the asymptotics of A L ((3,p > p m (/3)) 
is completely different than (11.511) and it is equal to 

A L (f3, p) = \p(p - pUPWl}" 1 + o(V^) . (1.53) 

Since Vl = L 3 e 2aL , we obtain: 



lim — 



L ^ V * M *x^i>M) ^m 2 *Vlm) J2U (k 2 - 1)/L] + A L ] - 1 
lim J- ln[/3A L (/3, p)]" 1 = 2e*/Aj = p m (/3) - p c (/3), (1.54) 

and there is a macroscopic occupation of the ground-state that yields for the term 
corresponding to k = (1, 1, 1): 

p - p m (/3) = L nm Y L e^(M-l > ° " (L55) 

Notice that for p c (/?) < p < Pmi.fi) we obtain the vdBLP-gBEC (of the type III), i.e. 
none of the single-particle states are macro scopically occupied, since the exponential 
anisotropy: L\ = L2 = Le aL , L3 = L, and (ll.5ip imply for any mode k: 

p k (/3, p) := Km Y L 6 ^UM)-1 = ° " (L56) 
On the other hand, the asymptotics A L ((3,p > p m {(3)) = [f3(p — p m (,S))Vy -1 implies 

P*(W)(P,p) == ^Y L e^~^M)-l = ° ' (L57) 

i.e. for p > p m (/3) there is a coexistence of the saturated type III gBEC, with the 
constant density (ll.54p . and the standard, i.e. the type I vdBLP generalised BEC in the 
single (ground) state (jl.55p . |BZ] . |MuSa] . 
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Since the inverse temperature (3 and the chemical potential p are the intrinsic param- 
eters of the Gibbs grand-canonical ensemble, it is natural to use them in our analysis of 
BRPF in the thermodynamic limit. However, as it is indicated above to ketch the BEC 
in this ensemble we have to make the chemical potential Independent. To elucidate 
this L-dependence for anisotropic boson systems we consider first the average density 
p(L, A), i.e., the Gibbs expectation value of the total particle number in reservoir 
divided by the volume Vj> 

p(L,A):=i^W(fc,L,A) , (1.58) 

L k6N 3 

where 

W(k, L, A) := = . (1.59) 

ex V [((3hHy2m)j: 3 j=l (k]-l)/q + (3A]-l 

Here and hereafter, we set > fixed, but arbitrary, and suppress it as argument if it 
will not produce any confusion. 

Recall that p = £i(L) — A and the condition A > is needed for definiteness 
(stability) of the ideal Bose-gas. Then for fixed A > ( or for fixed p < ei(L)), we 
obtain in the limit 

_, . _, 1 f dk 1 f dk _ _ 

P{ } := l™ P{ ' } = J^rf J R3 e^^ 2fc2 /2 m +A) _ i < (2^)3 J R3 ^y/ta _ i = pc 

and lim A;o p(A) = p c . 

To make p(A) > p c , we let A be L-dependent. In fact, if Hindoo A(L) = 0, we have 

li* J2 W(k,L,A{L)) = p e 

— irvi 1/ r f * 



L— >O0 Vj 

keN 3 ,fc 3 >l 



and 



lim -i- Yl W(k,L,A(L)) = lim -^-|w(l, L, A(L))+ ]T L, A(L))| 

kGN 3 ,fc 3 =l keN 3 ,(fci,fc 2 )^(l,l),fc3=l 

see (lA.ip . If we set A(L) = e _rL (0 < 7 ^ 2a), for example, we have 

p(A) = Jam p(L, A(L)) = — + p c < — + p c =: p m . 

In this L-dependence on A, we obtain the type-Ill gBEC, i.e. none of the single- 
particle states are macro scopically occupied, and there are macroscopic contributions 
to the averaged density from infinitely many but relatively small parts of single-particle 
states whose quantum numbers lay in (N 2 — {(1,1)}) x{l} out of N 3 . Note that the 
first term of (ll.60p is zero and the second term yields p(A) — p c . 
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If we set A(L) = (L 3 e 2aL 5)~ 1 (5 > 0), we get 

lim p(L,A(L)) =p m + 5, 

L—>oo 

which indicates that there is a coexistence of the saturated type III vdBLP generalised 
BEC and the standard (i.e. the type I vdBLP in the single-particle ground state, |BZj . 
[MuSap . Note that the first term of fll.60|) yields p(A) — p m and the second term yields 
Pm - Pc- 

(b) Spacial distribution of BRPF in anisotropic SLAB. Now let us examine the 
detailed position distribution of these three phases of the exponential SLAB system. 
The above observation indicates that it is convenient to introduce the parameters: 

"'^J^ gAWL^ - «.==«-^log[(L' e ")Aa)A(L)H (1.61) 

to distinguish these phases. We restrict our attention to those asymptotics for A(L), 
that give finite values of K\ and n 2 . Note that there are differences in factors in K\ 
and k 2 with those quantities in (11.601) . We use the present form to make the results in 
the following theorems simple. The differences in the factors occur because we consider 
the different object in the following argument from the above. It is the thermodynamic 
limit of the local density near the coordinate origin rather than the density obtained 
by averaging over the region A^. Correspondingly, there is a difference between critical 
values PcPm of the averaged density and p c , p m of our BRPF discussed below. In fact, 
we have p m ^ p m , while p c = p c . See Remark [1.91 and Remark [1.141 

First, we consider the microscopic distribution of the constituent bosons in terms of 
BRPF. The result is the following theorem: 

Theorem 1.6 // lini^oo A(L) = Aoo > 0, then the BRPFs { v^k l }l>o converge 
weakly to the BRPF z^a^o- If A^ = 0, then the BRPFs { ^,a £ }l>o converge weakly to 
the BRPF u 0tK for k = k\ + k 2 . Here ^a^o is the BRPF characterised by the generating 
functional: 

[ e-^u Aaofi (dO = Bet[l + Kf-]- 1 , 
Jq(r 3 ) 

with Kf°° = Vl - erf G(e^ Aoo -G)~ 1 Vl - erf. Whereas the BRPFu 0>K is characterised 
by the generating functional: 

[ = (1 + K{s/l-e-f, (1 + A / °)- 1 v / l-e-/))^ 1 Det [1 + Kj]~\ 

with K° f = Vl - e~f G{\ - GyWl - erf = Kf x=Q , where f G C (M 3 ) is non-negative 
and G = e^ h A / 2m is the heat semigroup on the whole space M 3 . 

Remark 1.7 In the case A ro = 0, the operator G(e /3A,X — G)^ 1 is unbounded. For a 
careful argument of this case, we refer to ITIblJ . 
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Remark 1.8 Let fx n , fj, K , jt4 be BRPFs whose generating Junctionals are given by 



[ e-tf>V(^) = Det[l + #°]-\ 
JQQS. 3 ) 

e~ m fiM) = (1 + e-/, (1 + ^j-Vl-e-/))" 1 

J QCM. 3 ) 



13) 

and 

■(/-9„H^1 = P ~r( yfl-e-f , (l+K° f )-i y/l-e-f ) 



TTien we /iai>e f/ie convolution 

Vq,k = P-n * Pk 

and the decomposition into ergodic components 

POO 

/ — t (ex) j, 

Pk = e p y Kt ; a!t, 

which shows that v n is a mixture of 'pure states p n * p^t (£ > 0). 27m mixture property 
would be grand canonical ensemble artifact. Then this is a manifestation of the Kac 
decomposition. 

Remark 1.9 From the above expression of the generating functional, we have 

(f,0 Vo,MQ = / Pf( x ) dx 



for p = K\ + k 2 + K Acc (x,x), where K A °°{x,y) is the kernel of the operator G(e^ A °° — 
G)- 1 , i.e., 

r e ik-(x-y) fifc 

K °°{x,y) = e/9( ^| fe | 2/2m+Aoo) _ 1 -ppja- (1-62) 

Note that K A °°(x,x) is x -independent. Hence, p may be regarded as the density of the 
system, which is uniform. This uniformity is the result of the thermodynamic limit. 

In fact, p corresponds to the local density of the compact region near the origin 
before thermodynamic limit, where the local density of the condensed part is expected to 
depend on the position because the ground state wave function of the Dirichlet boundary 
conditions is not uniform. When A^ > 0, K\ = k 2 = and p < p c . While p ^ p c 
corresponds to the A^ = case, where 

1 dk 



Pc = Pc 



e /3/l2|fc|2/2m _ l (27r) 



Remark 1.10 As one can see from the argument above and the proof of the theorem 
that K\ is a contribution of the ground state 0i,a, ^2 comes from the Riemann sum which 
is the sum of contributions of the states {0( s1iS2 ,i),a I s i, s 2 £ N, (si,s 2 ) ^ (1,1)} and 
K(x, x) is a contribution of other states. That is, if K\ ^ 0, the system enjoys a usual 
BEC (i.e., type I of gBEC). If k 2 ^0, it is in type III gBEC. WIT/ 
This feature also holds for Theorem and Theorem \1.15l 
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Remark 1.11 In the above argument, we have derived the "density" p from the view 
point that the system is specified by (3 and the L-dependence of A (i.e., p). However 
it may be natural to understand that the L-dependence of A is induced from and p. 
From this point of view, for (/3,p) in the region p ^ p c , we may take L-independent 
A = Aqo which is the unique solution of 



P 



j 



1 



dk 



g/^lfcp/am+Aoc) _ i (2tt) 3 ' 

for {/3, p) in the region p c < p ^ p m , A(p, L) is given by 



A(p,L) 



-f}-Kh 2 {p-p c {fi))L(l+o(l))/m 



where p m = 2ma/ (f3nh 2 ) ; 

for (f3,p) in the region p > p m , A(p,L) is given by 



A(p,L) 



8 



/3L 2 e 2 ^(p- Pm (/3)) 



;i +o(i)). 



Now let us consider the type of macroscopic distribution of the constituent particles 
in the exponential SLAB by using two different scaling arguments. 

First, we deal with the scaling, in which can be seen in the thermodynamic limit 
as a two-dimensional square S = [—1/2, 1/2] 2 , i.e. a finite infinitely thin slab. To this 
end we scale point measures by the transformation 



T S :Q(A L )3£ = J2 5 ^ 



V 



L 3 



a 2aL \x[ j) 



e M(S) 



where Ai(S) is the space of all (locally) finite non-negative Borel measures on 5*. 

Now let us introduce the following Random Field (RF), i.e., a Borel probability 
measure on A4(S) : 



(S) 
p,L 



-l 
S ' 



1.63) 



and also define the RF *M). h , which is characterised by the generating functional 



e- {f ' v) ^f(drj) = exp [ - a / f(x) dx 
M(S) L J s 



x 1 + b / f(x) ( Y\ cos 2 KXj) dx , 
Js j=1 

for non-negative / G C(S). Then we get the following result: 

Theorem 1.12 The random measure ^jfl converges to with a = K A °°(x,x) + 

k 2 /2 and b = Ki/2 weakly, with K Aco (x,x) defined by U.61\) and U.6ty) . 
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Remark 1.13 Let 5 g be the "deterministic" RF which has point mass on the singleton 
g dx G M.{S), i.e., 

[ e- {f ' v) 5 g (dri) = e~ & «*» . 

/m(S) 

Then, ^# a f) /ias i/ie decomposition to deterministic RFs: 

POO 

./o 

u>/iere 

2 

<7t(x) = a + bt( Y\ cos2 ^j) • 

3=1 

T7ie /irsi term o/ <ft a contribution of the mixture of the normal and the type-Ill 
gBEC phases, while the second comes from the usual BEC, see Sec \1.2\ Sec \1.3\ and in 
particular U.61\) . 

The Xi,X2 dependence of the second term reflects the square of the ground state wave 
function of the Dirichlet boundary conditions. In fact, it is the scaled and squared ground 
state wave-function after integrating the third variable x 3 . 

Remark 1.14 We have 

[ (f,v}^lf(dv)= [ P(x)f(x)dx, 

J M{S) JS 

where 

2 

p(x) = K + — + — || cos nxj 

3=1 

can be regarded as the density function of the system from the macroscopic point of view. 
Note that by averaging p, we have 

f p(x) dx/\S\ = K A -(x,x) + | + ^ = p, 
«/ s 

where \S\ = 1 is the Lebesgue measure of S. This, the averaging density of grand 
canonical ensemble U.60\) is re-derived. 

The second scaling we consider here concerns another macroscopic (mesoscopic) 
distribution, where the system looks in the limit as an infinite plane slab with a finite 
thickness: D = R 2 x [-1/2, 1/2]. 

Let us scale point measures by the following transformation: 

T D : Q(Al) 3 £ = ■— ► ^ = ^3 E W E M ^ C M W> ^ L64 ) 

3 3 
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where D L = [-e aL /2, e aL /2} 2 x [-1/2, 1/2] and M(D), or M(D L ), are the space of all 
locally finite non-negative measures on D, or D^, respectively. 
Let us introduce a RF on M.(D) by 



1.65) 



We also define a RF on M.(D) which characterized by the generating functional 



e-^^\dr)) = exp -a/ f(x)dx 
M(D) L 



1 + 6 / /(x) COS 2 7TX 3 <£c , 

where / G Cq(D) is non-negative. We get the following result: 
Theorem 1.15 TTie random measure ^(^L converges to 

a = K Aoc (x,x) and b = K\ + k 2 , 
weakly, with K(x, x) defined by (Q7 



Remark 1.16 Here k is same with that in Theorem \1.6l We can also get a similar 
decomposition of the RF as in Remark \ 1.131 

Remark 1.17 The "density" function of this scale is 

p(x) = K Aoc (x, x) + («i + K 2 ) COS 2 7T£ 3 . 

The second term is proportional to the value at (0,0, xs) of scaled squared ground state 
wave function. 

Remark 1.18 It would help understanding the whole picture of the SLAB system to 
summarize the results of the three scales roughly as 



-</,«> 



i/(df) = (1 + (p - p c ) + (y/l - e-f, (1 + Kf-rWl ~ e- f 



'Aool-l 



-^^ s \dri) 



xDet [1 + K l f 

(p A p m - p c ) + 



1.66) 



exp 



M(S) 



+ |0Ap c ) / f{x)dx 



1+ (P Pm)+ / f(x)(f[cos 2 7rx J )dxy\ (1.67) 

^ 5 7=1 



-^JZ( D \drj) 



exp 



M(D) 



(pAp c ) J f(x)dx 
l + (p-p c )) + / f(x) cos 2 n x 3 dx 
20 



D 



1.68) 



where c + = max{c, 0} and c A d = min{c, d}. The first equation expresses the BRPF 
for the microscopic scale, the second the BRPF for the macroscopic and third for the 
mesoscopic scale. Here, the parameters appeared in every scale are written in terms of 
p which appeared in microscopic version in Remark \l.*A 

From the mesoscopic scale ( fl. 68\) . we can see the £3- dependence of the local density, 
which is a effect of type-Ill gBEC. The term containing (p — p m ) + in the denominator 
in fli.ff7| ) shows the X2- dependence of the local density in this scale, which is a effect 
of the usual BEC. In this way, the qualitative difference between the two BEC phases is 
obvious. 



1.3 Microscopic and macroscopic distribution of constituent 
bosons in algebraically anisotropic Casimir prism (BEAM) 

(a) Casimir prism: one critical temperature and type-II gBEC. In this subsec- 
tion, we consider the distribution of constituent bosons in the Casimir prism (BEAM) 
|PuZ] in order to compare it to the SLAB case |BZ] . 

To review the picture of BEC in the BEAM system [BLP] or |Zae;2|, let us consider 



first the expectation value of the particle density in the grand canonical ensemble for 
the prism Al '■= L\ x L2 x L3, where L\ = L 7 , L2 = L3 = L. Then its volume is equal to 
Vl '■= \Al\ = L 2+1 . Similar to the case of SLAB, for a fixed A > 0, we get the limiting 
particle density 

p(A) := hm p(L, A) = JL ^ -^g^ _ { < p c , 

and limA^o p(A) = p c , where p c is the same critical density as in the SLAB case. 

To make p(A) > p c , we must let A be L-dependent: A = A(L). Note that in this 
case, we have 



hm-i Yl W(k,L,A(L))=p c 



keN 3 ,fc 2 >l,A:3>l 



if lim^oo A(L) = 0, see (1A.5j) . The remaining terms are 



y Yl W{k,L,A{L)) = ±.{w{l,L,A(I<))+ E W(k, L, A(L))| 

L keN 3 ,fc 2 =fc3=l L keN 3 ,fci>l,fc 2 =fc 3 =l 

-i^km + ^ 0{LlAA{Lrl,2h (L69) 

see (lA.llj) . where we put L 1 instead of L. Here the first term in (11.691) is a contribution 



from the single-particle ground state, and the second term is a contribution of countably 
many single-particle excited states. 

Notice that for 7 < 2 the lim^oo p(A(L)) =: p > p c implies A(L) = 0(L~ 2-7 ) = 
0(1/Vl). Hence, only the first term in (11.691) is non-zero in the limit, that yields the 
type-I (ground-state) BEC. 
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To insure that lim^oo p(A(L)) =: p > p c for 7 = 2, one must put A(L) = 0(L~ A ) = 
0(1 /Vl). Then (I1.69P implies that not only ground-state but also the contribution to the 
second term from each single-particle state is nonzero in thermodynamic limit. Hence, 
for 7 = 2 the gBEC in BEAM is of the type-II [BLP] . [Zag2]. 

If 7 > 2, then to insure lim^oo p(A(L)) =: p > p c we have to put A(L) = 0(L~ 4 ) = 
0(l/y L 4/(2+7) ) [BLP] . It is known that in this case there is no macroscopic occupation 
of any single-particle state, which corresponds to the type-Ill gBEC, [BLP], Zag2 . 



In the following we consider only the most interesting boundary case 7 = 2, i.e. the 
case of the type-II gBEC. 

First, we examine the microscopic distribution of constituent particles in this kind 
of the Casimir prism (BEAM) in the thermodynamic limit. We consider also z^,a £ and 
the BRPF (HTTP]) for A = A L and in the limit L -> oo. 

We introduce a function <p on [0, oo) defined by 

oo 

: = E (2B -!).-! + „ ■ < L7 °) 

n=l v ' 



It is not difficult to show that (p is decreasing and its asymptotic behavior is 

(fix) 



l/x + 0(l) for small x, 

7t/(4a/x) + 0(1 /x) for large x, 



(see Remark 1 1.25ft . For the present system, we use the quantity 

«=-^limv,(2mL 4 A(L)/7r^ 2 ). 

We will see that each term of </? reflects the contributions of each state in { (f)( Slj i t i)A | si : 
odd}, in the proof of the theorem. This indicates that the type II gBEC occurs in the 
system for the case k^O. 

Now we can formulate the following theorem: 

Theorem 1.19 If lim^oo A(L) = Aoo > 0, then the BRPFs { ^,,a l }l>o converge 
weakly to the BRPF z/ Aoo ,o- If lim^oo A(L) = 0, then the BRPFs { 

V[_i : a l }l>o converge 

weakly to the BRPF Vq^. Here the BRPFs : va,o, ^o,k (with k replaced by k), coincide 
with those appeared in Theorem \1.6l 

Remark 1.20 As in Remark \1.9[ we have 

p = R + K Aoc (x,x) 

for the expected "density" of the system. 

We may also derive L-dependence of A(L) in order to give a prescribed (3 and p. 
That is, 

for (f3,p) in the region p ^ p c (normal phase), we may take L-independent A which is 
the unique solution of 

1 dp 



P 



; /3(7> 2 |p| 2 /2m+A) _ 1 ( 27r )3 
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and for (/3,p) in the region p > p c (type-II gBEC case), A(L) is given by 

n 2 h 2 „ l( {p-p c ){3ir 2 h 2 ' 



y ' 2mL 4 ^ V 16 / 



2mL^ V 16 
where p c is the same one as in the SLAB case. 

Now let us consider the two type of macroscopic distribution of the constituent 
particles of the Casimir prism by using two different scaling arguments. 

In the first scaling, looks like i? = Rx [—1/2, 1/2] 2 and the corresponding trans- 
formation is 

T R : Q(Al) 3 £ = 5>. w — ► v = ys e M(R). 

j j 

Let us introduce the random field 

= "Mi (T R )-\ (1.71) 

( Fi) 

We also define RF ^# o 6 , which is characterised by the generating functional 

e- {f ' v) ^f(dr]) = exp [ - a [ f(x) dx 

m{r) l y B 



1 + / f(x) ( cos 2 7rXj) dx 

J R „■ n 



i=2 

for non-negative / G Cq(R). We get the following result: 
Theorem 1.21 The random measure converges to ^ ah with 

a = K Aoc (x 1 x) 

and b = k, weakly. 

Remark 1.22 This distribution is a convex combination of "deterministic states" as in 
Remark[TTM 

Remark 1.23 We have 

(f,Tl)J$(dT,)= f P(x)f(x)dx, 
M(R) JR 

where 

3 

p(x) = K Ao ° (x, x) + k Y[ cos 2 nxj 

3=2 

is the "density" of the system in the mesoscopic scale. The second term of p is pro- 
portional to the value at (0, X2,x%) of the scaled and squared low-lying eigenfunctions 

{ 0(s,O,O)A L KeN- 
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In the second scaling for Casimir prism, the system looks like I — [—1/2, 1/2] in the 
limit. The transformation for the scaling is 

Ti : Q(Al) 3 t = J^u) ,— ► T] = jjiYtty/ifi E 
j j 



We introduce the RF 

We also define the RF j&V> on A4(I), which is characterized by the generating func 



tional: 

exp - if A °° (x, x) J T f(u) du 

e-^^\d V ) = 1 ntM4 _ fp1 1 , (1-73) 

'M(i) Det [1 + /it] 

where R is the integral operator whose kernel is 

, 8m cosn7r(w — i>) — cosrarfl + u + v) , „ . 

*«■«)== ^ E — 5 — 1 (174) 

n=l 

8m 7r [cosh (7ra(l — \u — v\)j — cosh (ixa{u + t>))] 
flh 2 7r 2 2a sinh 7ra 

with 



;2mL 4 A(L) 
a = r hm V X ' 

L->oo V fi - i 7T z 

and / G C(J) is non-negative. The second equality of (jl.74p is derived by the Fourier 
series expansion for cosh. Note that a may be pure imaginary. We understand that 
R = for a = oo. We get the following result: 

Theorem 1.24 The random measure converges to weakly. 

Remark 1.25 We have 

(f, V )^ I \d V )= [p(u)f(u)du, 
M{I) Jl 



where 



. . a . . „, , a . . _ cosh(7ra) — cosh(27ran) 
p{u) = K A °°(x,x) + R(u,u) = K A °°(x,x) + k- v ; 



cosh(7ra) — 1 

is the expected "density". Note that (p(a 2 + 1) = R(0,0), R(0,0) = n and 

16m 



7T COsh TTy/x — 1 — 1 

<p[x) = — -j= . — 

Ay/x— 1 sinh7rvx — 1 

holds. 
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Remark 1.26 In this case, the macroscopic "density" has random distribution. In fact, 
it follows from |T773|) that the density is a squared Gaussian RF. The operator R consists 
of the contributions of low-lying eigenfunctions { 4>( s ,o,o)A L KeN- This revived random, 
which is disappeared once in mesoscopic scale, may be considered as the contribution 
of boundary effects through the law-lying eigenfunctions. This are strong fluctuations, 
which are typical in type-II gBEC, see e.g. jBLLf . JPuZf . 



2 Proofs for the SLAB 
2.1 Proof of Theorem [1J3] 

The function /, which will appear in the generating functional Z\(f) is assumed to be 
non- negative continuous function on IR 3 with compact support. In the course of the 
thermodynamic limit, / is regarded as a function in Cq(Al) naturally. 

Let P be the orthogonal projection operator on L 2 (A L ) onto the closed subspace 
spanned by {</>( Sl)S2 ,i)A I s i> s 2 £ N}, and Pi the orthogonal projection onto j£f = J^ 1 - 
Let us recall that the function K\(x, y) is given by (jl.2ip with L\ = L 2 = Le aL , L 3 = L 
and L-dependent A = A(L), though the dependence on A(-) is suppressed on the 
symbol. In the following, we also suppress the Aoo-dependence from K A °°(x, x), 
etc., and we write A instead of A(L). 

Lemma 2.1 The following convergence of the operators holds in trace class: 

^/l-e-fPJtAP^l-e-f — ► Kf 



Proof : From ( 1A.4I) , 



oo oo 



P 1 K A P 1 (x,y)= E 



4>s,A(x)4> s ,A{y) 



si,s2=i S3=2 exp 



2m I L 2 e 2aL ^ L 2 



f 

Jm 



e 2nip-(x-y) 



exp 



dp = K(x, y) 



holds locally uniformly in x, y. Since supp / is compact, we have 

Tr [y/l - ~^P X K K P 1 \T\ - e~f] = [ P^aP^x, x)(1 - e~ f{x) ) dx 

J A 



On the other hand, 



■» / K(x,x){l-e- f(x) )dx = TrK f . 
\fl - e-f PiK K Pi \A ~ e-f -K f \\ 
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= sup (y/l-e-f<p, (PiKaPi - G(e@ A - G) -1 ) \A - e~f<ip) 

V ,ip£L 2 {A) 

^ sup / y/l - e-fW\tp(x) | y/l - e~fW \^{y) \ dxdy x sup \P 1 K A P 1 (x,y)-K(x, 

tprt J Ax A £,3/Ssupp/ 



^ sup || 1 — e" 



-/(*) I 



1IMI2IMI2 sup \P 1 K A P 1 {x,y)-K{x, 

z.ySsupp/ 







holds. Then, the Griimm convergence theorem Zag| yields the lemma. □ 
Lemma 2.2 TTie following convergence of the operators holds in the trace-class topology: 

P K A P — ► K y/l-e-f(y/l-e-f, •) . 
Here the limit operator ny/l — e~f(y/l — e~f, ■) is a projector, which acts as 

ip 1 — > (ki + k 2 )(v / 1 - e _/ ,V>) \A - e~ / , « = «i + k 2 , 

/or^ e L 2 (M 3 ). 

Proof : In the expression 



P ^aPo(x,2/)= ^ 



0(-i,M,i) 1 A(a?)0( w ,«i,i) ) A(3/) 



si,s2=i exp 
the numerator may be represented as 



2m L 2 e 2aL ^ 



0( Si , S2 A),a(z)</>( S1)S2 ,1),a(Z/) 



(2.i; 



n 



Le aL L 



— — h (-1) 3 sm 2 J „V / J -sin 2 n 3 yv 



3=1 

Hence we get 



2Le aL 



2Le 



aL 



\ 2 ttx 3 7ry 3 
> — cos — — cos — — . 

J L Lj ±j 



P K A P (x,y) = Yl 



cos(7ix 3 /L) cos(7ry 3 /L) 



JJ, e 2aL 
si,S2SN,odd G X P 



where 



1^1 ^ JJ, e 2aL 



2m L 2 e 2aL ^ 



C(s 2 + s 2 )/L 2 e 2QL 



exp 



r 



where C depends only on \x\ and \y\. It follows from (1A.2I) and ( 1 A. 81) that 



P ^APo(x,t/) 



/3AL 3 e 2aL (3nh 2 L 
fO(L- 1 ) + 0(L- B e- 4aI '( / 9A)- 1 



(2.2) 



Ki + K 2 



(2.3) 
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locally uniformly in x,y when L — > oo. The first term of the 2nd member in the above 
equation corresponds to s = (1,1). Now, the lemma can be proved by the argument 
similar to the proof of the previous lemma. 

□ 

By consequence the weak convergence of the BRPFs is implied by the following 
convergence of the generating functional: 

/ e-^v^dg) = Det [1 + v 7 ! - e-f(P K A P + Pi# A Pi) V 7 ! - e~f ] 1 



-i 



— > Det [1 + «V 1 - e~f{ V 1 - e-f, +K f ] 
= Det [1 + K^Det [l + Ky/l - e~f( y/l - e~f, (1 + Kf)' 1 ]' 1 
= Det [1 + K f ]- l {l + K(y/l-e~f, (1 + K f )~ Vl - e~f ) } _1 

= [ e-^uM) ■ 

Here we used Lemma 12.14 Lemma 12.21 and 

Det [1 + A + B] = Det [1 + B] Det [1 + A(l + By 1 ] . 

One has also to notice that projector: k\/1 — e~f (a/1 — e~f, •) , is one-dimensional op- 
erator. 

2.2 Proof of Theorem [LT2] 

Let P be the orthogonal projection operator on L 2 (Al) onto its one-dimensional sub- 
space C0ia l , and Q the orthogonal projection onto (C0ia l ) . 
Put 

f L {x) = L^e- 2aL f{x l /Le aL ,x 2 /Le aL ) 

and 

f ( - L \x) = L 3 e 2aL (l-e~ f ^ L3e2aL ) 

for bounded measurable non-negative function / on S — [—1/2, 1/2] 2 . The functions 
/l,/^ are defined on and on S, respectively. Note that f^ L \x) t f(x). Then we 
have 

(h,0 = (f,v) = (f,TsO- (2-4) 
From ( ITTOj) . pg| . (E2D and Q3| . we have 

/ e-^^ S l(d V ) = [ e-^v^M) = Det [1 + K\ L (1 - e^)]" 1 . 

JM{S) JQ(A) 

By virtue of the Feshbach formula we get 

Det [1 + K A (1 - e"^)]" 1 = Det [1 + (P + Q)tf A (l - e"^)(P + Q) ]~ l 
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Det [l+PK A (l-e- fL )P-PK A (l-e' fL )Q(l+QK A (l-e- fL )Qy 1 QK A {l-e- fL )P]' 1 

xDet[l + QK A (l-e- fL )Q}-\ (2.5) 

Note that 

PK A P = (e pA - 1) _1 P, \\PK A P\\ < ((3A)-\ 



\\QK A Q\\ = exp 



(3H 2 ti 2 {A-1) 



+ /3A 



l) =0(L 2 e 2aL A^A)- 1 ) 



and 

It follows that 



2mL 2 e 2aL 

l-e-^|| ^ \\f\UL s e 2aL , \\QK A {l-e-^)Q\\=0{L- 1 ). 

\PK A {1 - e~ fl -)Q(l + QK A (1 - e- fL )QY l QK A {\ - e~ fL )P\\ 
1 



: o(l- 1 ) = o(l- 1 ; 



/3A L 3 e 2aL 

Here, we recall the assumption of boundedness of (AL 3 e 2aL ) _1 and that QK A (l — e~^ L )Q 
is positive. The convergence of the one-dimensional operator 

PK A (1 - e-f»)P = (^A,a-^)^) p 

2 



eP A - l)L 3 e 2aL Js 
1 



f {L \x)( JJcos 2 TTXj) dxP 



lim - 

l^oo f3AL i e 2aL 



3=1 

2 



/(x) ( J^J cos 2 7TXj) dx P 



3=1 



holds. Note that the trace norm is identical to the operator norm for a one-dimensional 
operator. Thus we get 



Det [1 + PK A (l - e~ fL )P- PK A (1 - e~ fL )Q(l + QK A (l - e' fL )Q) QK A (1 - e~^)P\ 

4 



lim - 

l^oo (3AL 3 e 2aL 



f(x) ( TT cos 2 -KXj) dx. 



(2.6) 



3=1 



We remark here that the function Ylj=i cos2 7lx j ls a reminiscent of the ground-state 
wave function. 

Next, let us consider the operator QK A (1 — e~f)Q. To simplify the notations, we 
introduce the following symbols: 



Sl,S2 



/ (L) (M)[n cos ( 27rs i M i)] du 

3=1 



(2.7) 
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the corresponding limit: 



f(u) [ J cos(2irsjUj)~\ du 

3=1 



(21 



for s = (si, s 2 ) G N 2 , where N = {0} U N, and the symbol: 

-I3h 2 n 2 / S 2_i + S 2_i s 2_ x 



W(s, L) = exp 



. 2m 



L 2 e 2aL 



V- 



/3A 



- 1 



-i 



(2.9) 



for s = (si, s 2 , s 3 ) G N 3 . 

Then argument based on the Riemann-Lebesgue lemma yields the following result: 



Lemma 2.3 For every s = (s±, s 2 ) G N 2 , 



K-ai L) l < 



/2L 3 e 



3„2aL 



|ai L )|<6. 



holds. And a* converges to 0, i.e., 

Ve > 0, 3N G N, 3L > : si + s 2 > N , L > L = 
Returning back to the proof of Theorem 11.121 note that in the expression 

Ti[QK A (l-e' fL )Q} = Y,W(s,L) [ |0 S;A (x)| 2 (l - e"^) dx 



e / / w («)n(i - ™ 2 -w) *• 

seN 3 J5 j = l 



10,0) ^3 e 2aL ^ ^ 



Ri 



sen- 5 



eqs. dAT]) . (jO)l yield 



E -+E 



lim 



S^l,3 3 = l 



7T 



A(3hH 2 /2mL 



1 

+- 



log (L 2 e 2aL A (/3A)- 1 ) 



8 7 R 3 exp [/m 2 7r 2 |p| 2 /2m + /3A (£=oo) ] - 1 

as L — y oo, and we shall show that lim^oo i?^ = in Lemma [2.41 below. 

Then together with observation: ||Q.Ka(1 — e~f L )Q\\ as L — > oo, we get 

Det [1 + QK A (1 - e~^)Q] = e ^[QK A (i-e-fL )Q]+0 (i) 
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exp 



/ 777 

o ( o,o) (K(x, x) + — — lim - log (L 2 e 2aL A (/3A)- 1 



l32nh 2 l^oo L 

= e aJ s f(u)du (2.10) 

As a consequence we obtain from (12.51) . (12.61) and (12.101) that 

/ e-Vrt J(fl{&r]) e a Is f ( u)du (l + b [ f (u)(TT cos ttu A duY 1 
Jm(s) ' v Js jJi ' 



'M(S) 

This proves the Theorem 11.121 □ 
Now we return to the promised Lemma 12.41 



Lemma 2.4 



lim R L = 



Proof : With help of (12.71) and (12.91) we can re-write Rl as the sum of three terms 
R L = rM + Rf + Rf\ where 

V (-l)^ +1 o L W ^ S ' L) R {2) - V ( 1^+V W ^ S ' L) 



i ~ Z^ ^ ^ "("1.0) ^3 e 2«L ' Jt £ _ Z^ ^ L > "(0' s 2) £3g2c*L 

and 

For an arbitrary but fixed e > 0, let Lq, No be the numbers which appeared in Lemma 
12.31 Suppose L > L . Then, we have 



l/?«l< V \n L \ W{ - S ' L > 



< I O(l,0) | 



W(s,L) 

l"(si,0)l 

sGN 3 ,s^l 

Z^ L 3 e 2aL l a ( s i>°)l Z^ L 3 e 2aL 



L 3 e 2cti 

si=JVo+l s 2 ,s 3 6N 



+ E l a fi,o)l I] 



^[0(e^)+0(logL)]+(AW)||/IU^ 



°°Le aL 

W(s,L) 



p 2aL 

sSN 3 ,s^l 



O(A^ logL/L) + O(e). 
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We have used (1A.7|) at the third inequality and ( 1A.6j) at the last line. 

Hence by taking L large enough, we can make arbitrarily small. Thus we get 



lim^oo Rf = 0. 



We get limi^oo R K L J = in the same way. For \imL^, OQ R K L J = 0, we proceed in a 
similar manner, i.e., 

ir( 3 )i< \n L i ^ (S;L) 

\ L I ^ 2^ l a (*i,«a)l £3 e 2aL 



,^ H/((l,l,s 3 ),L) ^ ^ L l W((s u s 2 ,s 3 ),L) 

^l a (l,l)IZ^ L 3 e 2aL ^ Z^l^i^)' L 3 e 2aL 

s 3 =2 3<si+s 2 <iVo 83=1 



W(s,L) 

JJ, e 2aL 

S1+S2>Nq S3 = l 



oo oo 

Yl i a f»i,«)iz^ 



^ W((l,l,s 3 ),L) f N (N -l) x f H^UU) 



oo 

S3=2 S3=l 



= O(iV 2 /L) + O(e), 

where we have used ( 12. 8p and ( lA.lip in the last line. □ 
2.3 Proof of TheoremEE] 

Let L> L = [-e aL /2,e aL /2] 2 x [-1/2, 1/2] and D = M 2 x [-1/2, 1/2]. Put 

/ L (x)=Zr 3 /(x/L) 

and 

/( L )(x) = L 3 (l-e"-«^ L3 ) 

for non- negative / G Cq(D). For large L > we can naturally consider / G Cq(Di). The 
functions /l,/^ are defined on and on respectively. Note that f^ L \x) /(x). 
Then by definition (jl.64p we have 

(f L ,0 = {f,T D . (2.11) 

Recall that Po be the orthogonal projection operator on L 2 (A^) onto the closed 
subspace spanned by {(f>( Sl ,s 2 ,i)A I s i> s 2 £ N}, and Qi the orthogonal projection onto 

From ([Em (fl~T8D . (EH]) and (TL65D . we deduce that 

JM{D) JQ(A) 



31 



1-1 



= Det [1 + \/l - e~^K AL \/l - 
Det [1 + y/l-e-^(P Q K AL P + P x K Al P x )^/\ - e^]' 1 
= Det [1 + x/l - e~^P l K AL P 1 y / l - e^]" 1 



xDet [1 + y/l - e-^P K AL P ^l - e -^(l + v 7 ! - e-^PJi^P^l - e^)" 1 ]" 1 . 
We skip the proof of the following obvious operator-norm estimate: 
Lemma 2.5 

|| y/l - e-hP x K KL P x y/\ - e-/i|| = 0(L _1 ). 
On the other hand we need the proof of the following limit: 
Lemma 2.6 

Tr [y/l - ~/F^P x K Al P x \T\ - e~h\ > K(x,x) [ f(x)dx (2.12) 

L^oo J D 

Proof : By virtue of 

2 sin 2 + /-//,) = 1 - (-1)" cos {-j^Lyj), 
we get the following representation for the left hand side of H2. 12[) : 

3 



eN 3 iS3 ^i exp ^j=i + P A \ -l JD j=i J 



E 



J D f m (y)dy 



T3 P 2aL 

seN 3 ,s 3 ^i exp 
where Li = L 2 = Le aL , L 3 = L and 



2m 2-<j=l L 2 



/3A 



(2.13) 



y — 



(2.14) 



// (L) (y)[-E(-l)^os(^ V pL yj )+ J2 (-lY j+Sh cos (^%) cos {^Ly k 

-fK-lJ^COS (^/-// ; )_ dy. 



Now, let us put 



9l(p) = \ c 




2m ^J=l ^2 rpii 



for p 3 e (0, £ 



j = 1,2,3, S 6N 3 )S3 ^1 
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and 



a L (p)= I / (L) (y)[-E(- 1 ) Sj C0S (^^)+ E (-IJ^oos^Ly,) cos (^Ly 

" D i=l J l<j<fc^3 J J 



JJ(-l)^ cos (^^Ly^j dy for Pj - = ( 



Then we have 



1 



j = 1,2,3, sef. 



because 



3 s 2 - 1 3 s 2 
E^^Eo/a fOT '^,^1. 
3=1 ' 3=i L i 



We also have lim^oo a^ip) = a.e. by the Riemann-Lebesgue lemma, and 

lim ^I(P) = OH 91 19 /n ^.C. 

L-j-oo e /3ft 2 7r 2 |p| 2 /2m+/3A _ J 

Therefore, the dominated convergence theorem yields for (12.141) 



lim tl = lim / aL(p)^ l(p) dp = , 
and for the first term of the left hand side of eq. (I2.13P 

.L flL)(v) iy Lr * l(p) * -> L m dy L> e^Z^-i ■ 

Remark 2.7 From the above two lemmata we obtain that 

lim Det [1 + y/l - e-KP x K KL P x y/\ - e~fA = e ^,x)J D f( y )d y . 

We refer to hi. 62$ for definition of K(x,y). 
Lemma 2.8 

lim Det [1 + y/l - e-hP K AL P Q y/l - e~h{l + \A - ~e Z ^P x K Kh P x ^f\ - 



□ 



1 + K \ f( x ) c °s 2 7TX 3 dx 
Jd 
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Proof : Let Ul : L 2 (A L ) — > L 2 (D L ) be the unitary operator given by {Ui,g){x) 
L 3 ' 2 g{Lx) for x G D^. Then, the transformation law of integral kernels is given by 

(U L KU£ 1 )(x,y) = L 3 K(Lx,Ly) 

for x,y G Dl. Together with eq.( 12.1l) . it follows that 

L- 3 (U L P K AL P U^)(x,y) = P K A P (Lx,Ly) 

8L~ 3 e~ 2aL cos(7rx 3 ) cos(iry 3 ) 



si,S2SN,odd ex P 



2m L 2 e 2aL 



+ /3A 



Ql 



(2.15) 



with 



I ,< V- g (^ 2 (|x| 2 + |y| 2 )( gl + ^)/L 2 e 2 ^) 1/4 



si,s 2 sn 



exp 



2m L2i2SI "I" 



N 2 +l2/| 2 )0(L- 



-l/2\ 



where C is a numerical constant. We have used sin 2 a: «C C'lxj 1 / 2 and (1A.9j) . By (1A.2[) . 
we get the convergence 

L~ 3 (U L P K AL P U£ 1 )(x,y) — ► k cos ttx 3 cos 7ry 3 



uniformly in x, y on compacta. Then, by the same argument as in Lemma [2.11 we get 



U L Vl - e-^P K A P y/l - e-f^UZ 1 = y^fW ' L~ 3 U L P Q K Al P UI x 



— > Ky /cos7rx 3 (y fC0S7TX 3 , 

in trace class. Then the lemma follows from Lemma 12.51 and the fact that a unitary 
transformation conserves the Fredholm determinant. □ 
Remark 12.71 and Lemma 12.81 with arguments developed at the beginning of this 
subsection, imply the proof of Theorem l 1.151 



3 Proofs for the Casimir prism (BEAM) 
3.1 Proof of Theorem EES 

Recall that below the function / G C*o(lR 3 ) is non-negative and we consider / G Cq(Al) 
for restriction Aj, C I 3 . 

Let Qo be the orthogonal projection operator on L 2 (A L ) onto its close subspace J£o 
spanned by {0( Si ,i,i)a l | si G N}, and Q\ the orthogonal projection onto J%\ = (j^o) ± . 

Lemma 3.1 The following limit: 

Vl-e-fQ^AQ^l-e-f — > K f . 
is true in the trace-class topology \\ ■ 
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The proof of the lemma follows the lines reasoning of Lemma 12.14 using (1A.5[) instead 
of (TjOD. 

Lemma 3.2 One has the limit: 

|| ■ ||i - lirnQo-KAQo = «yi - e~/(v / l - e~f, , 
where k is defined by K3.2t) . 

Proof : Similar to eq. (I2.2p in the proof of Lemma I2.2[ we get 



n is n i \ ST 8 cos(irx 2 / L) cos(iry 2 / L) cos(irx 3 / L) cos(iry 3 / L) g 
Q K A Q (x,y) =: > — = hi<L , (3.1 

, l6 ^dd L exp ^£izi + /3A -1 



where 



\r l \ < £ 



7"4 

SlG N exp 



2m L 4 



0(L- 7 (f3k)- l ) + 0(L- 2 ) 



where C depends only on \x\ and \y\. We used flA.lOj) in the last estimates. Then, by 
virtue of ( 1A.3[) . we get that 



Q K A Q (x,y) 



16m 



lim (f(2mL i A/7r 2 h 2 ) = k 



(3.2) 



converges uniformly in x, y on supp /. The rest of the proof is similar to that of Lemma 

EU □ 

Now, taking into account Lemma 13.11 and Lemma I3.2[ the Theorem 11.191 can be 
proved along the same lines of reasoning as Theorem II .61 



3.2 Proof of Theorem[L2l 

Let R:=Rx [—1/2, 1/2] 2 and R L := [-L/2, L/2] x [-1/2, 1/2] 2 . Put 

f L (x) = L~ 3 f(x/L) , 

and 

f L \x) = L i {l-e-^l L3 ) , 

for bounded measurable non- negative function / on R. Then, the functions /l, are 
defined on and on R, respectively. Note that f( L >(x) f f( x ) f° r L — > oo. Then we 
obtain 

(h,0 = (fiV) = (f,T R , (3.3) 

see fll.64p . From (11.101) . (I1.18p . (13.31) and (ll.7ip . similar to arguments we used above it 
follows that 

Jm(r) Jq(a.) 
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= Det [1 + v 7 ! - e~^QiK kL Q x \/\ - e~^] 1 

xDet [1 + yfl - e-fiQ K AL QoVl ~ e~^(l + \A - e^QiK^y/l - e^)" 1 ]- 1 . 
Again as above, the following estimate in the operator norm is obvious: 

Lemma 3.3 

yi-e-^QrK^Q^l-e-^W = 0(L~ l ). 
We also have the following limit: 
Lemma 3.4 

Tr [y/l - e-^Q x K Kh Q x ^\ - e-h] ► K(x, x) [ f(x) dx 

L^oo J R 

Proof : Similar to proof of Lemma 12. 6\ we get 



y 1 



hf (L) (y)dy 



L 4 

s6N3,(si,s 3 )^(i,i) exp 



2 m 



where 



R 



L ■'- 



y - 



3 

1 



Rt 



(3.4) 



' ' L 4 \ /3h 2 -rr 2 , n A 

sen 3 ,(s2,s 3 )^(i,i) exp | }_^ j=1 -^3- + p A 

2tts 



x 



// (L) (y)[-E(- 1 ) S3cos (^P^)+ E 

-fi(-i) 



l<j<fc<3 
'27TS., 



ifc cos ( —j-^-Lyj ) cos 



27TS fc 



dy. 



In this case, we put 



cxp 





2m 2^3=1 ^2 



for pj G 



forp 2 ,Ps e (0, 



j = 1,2,3, sGN 3 , (53,83)^(1,1) 



and 



(P):= / / (L) (l/)[-E(- 1 ) Sj cos (^P^)+ E 

R .7=1 J l=£i<fcsS3 



ri(- i ) sj 



2irSj 



dy for 



27TS, 



ifc cos ( ^^-Lyj ) cos 



2vrs fc 



j = 1,2,3, sef. 
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as before. Then we have 



and lim^ooa^ (p) = 0. Therefore by the dominated convergence theorem, we get 

lim R L = lim / af \p)^f\p) dp = , 



'(0,oo) 

and for the first term of the left hand side of eq. (13.41) 

dp 



fV(u) du / ¥ L R \p) dp / /(«) / -^^fc^ • (3.5) 



□ 



Remark 3.5 From the above two lemmata we obtain 



lim Det [1 + y/l-e-f^K^QWl-e-h] = e K i*.*)J R m*> 



Lemma 3.6 



lim Det [1 + y/l - e-^Q K AL Q y/l - e~^[l + \A - e^Q^^Q^l - e^)" 1 ] 

L— >oo 

Det [1 + bsffl ® COS 7TX 2 (C0S 7TX 2 , ® cos 7rx 3 (cos 7rx 3 , a/7 •))] 

1 + b / /(?/) COS 2 7T?/2 COS 2 7T?/3 dy . 



Here, the Fredholm determinant in the first member is for operators on L 2 (Al), the 
second determinant is for operator on L 2 (R). Here I is the identity operator on L 2 (M.), 
which is considered as a component of the space L 2 (R) = L 2 (M) tg) L 2 (— 1/2, 1/2]) ® 
L 2 (-l/2,l/2]). 

Proof: Let U^' : L 2 (A L ) — > L 2 (R L ) be unitary operator defined by (Ui,g)(x) := 
L 3 / 2 g(Lx) for x G Rl-,9 £ L 2 (A^). Then, the corresponding integral kernels trans- 
form as 

{U^KijJ^y^y) = L 3 K(Lx, Ly) 
for x,y £ Rl. Then it follows that by (1A.12j) 

L-^U^QoK^QoiU^y^v) = QoK A Q (Lx,Ly) 
cos(iiX2) cos(7iy 2 ) cos(nx 3 ) cos(7iy 3 ) . nsi nsi . .irsi irsi 



2-77 r ^ ] sm(— +— L*)sm( — +— Lx) (3.6) 



7"4 

Sl =i exp 
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oc 



8m v-^ COSO — C0S7TS! 

2^ 2 , — cos (™ 2 ) cos(7ry 2 ) cos(7rx 3 ) cos(7ry 3 ;, 



= K COs(7r£ 2 ) COs(7n/ 2 ) COs(7nC 3 ) COs(7n/ 3 ) 

locally uniformly in x, y. Recall that here we used 

16m /2mL 4 A|L_s.oo' 



K :- 



and definition (jl.70p . 

Then, we obtain that the II • Ih-norm limit 



— y ka/// ® cos(7TX2)(cos(7TX2), ® cos(7rx 3 )(cos(7ra;3), a// •)) ; 

exists. Thus, we get the lemma. □ 
Remark 13.51 and Lemma 13.61 yield Theorem ll.21[ 

3.3 Proof of Theorem[L2! 

For non-negative continuous function / on I = [—1/2, 1/2], we define 

f L \u) :=L 4 (l-e-^/ L4 ) , 

and 

f L (x) :=L- A f( Xl /L 2 ) , 

for x = (xi,X2,Xs)- The functions /l,/^ are defined on and on J, respectively. 
Note that / (L) t / holds point -wise for L — y oo. 
By ffTMj) we have 

(fL,Z) = (LT I 0, (3.7) 

and (ITTUj) . ffTTT8"D and (jT772]l yield 

M{I) JQ(A) 
= Det [1 + y/l - ~^QxK Kl Q x sT\ - e~^]~ l 

xDet [1 + \A ~ e-fLQ K AL Q Vl ~ e~^(l + \Tl - e I hQ 1 K KL Q x <s/\ - e-/*)- 1 ]- 1 , 

similar to what we established above. 

Again the following estimate in operator norm is obvious: 

Lemma 3.7 



yi - e-fiQ 1 K AL Q 1 Vl - e-ft\\ = 0(L~ 2 ). 
We also have the following limit: 



38 



Lemma 3.8 



Tr [y/l - e-hQ x K hL Q x >/\ - e 1 ^} > K(x,x) [ f(u) du 

L^oo Jj 

Proof: Similar to the proof of Lemma 12. 6[ we get 

Tr [y/l-e-fLQ^&y/l-e-fL] 
1 fjf^H^du 



E 



L 4 



seN3,( S2 , S 3)^(i,i) exp 



/3h 2 n 2 ST 3 i _1 



E =1 -V + ^ 



2m £-^3=1 L 



Rj 



(3.8) 



where 



R 



L ■'- 



E 



L 4 



X 



ex P I Ej=i ~l| 

// (L) (l/)[-E(- 1 ) S3cos (^^)+ E (-1)^008 (^24/,) cosily, 

j=i i 

In this case, we put 







2m Z-ij=l L -2 



-0A 



for pj G 



j = 1,2,3, seN 3 , ( S2 ,s 3 )^(l,l) 



forp 2 ,P3 e fo, i 



and 



W = // (i) (y)[-E(-l)^os(^^)+ £ (-1)"^ cos (^Ly.) cos (^Ly k ) 



-IK- 1 )' 

Then we have 



cos 



i=i 



27TS 7 - r 



dy for pj- e ( Sj l - , 



j = 1,2,3, seN 3 . 



^(p)K$(p) = e ^ 2H ; /4m _ i , $ GL i(( ,oo) 3 ), 



km * L (p) e ^ ft 2 7r 2 b | 2/2m+ ^ A _ ! 



a.e.. 



and limx^jx, (p) = 0. The dominated convergence theorem yields 



lim 2?£ = lim 

L—>oo L— >oo 



(0,oo) 3 
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«l(p)^l(p) dp = 



and 

dp 



f {L) (u) du / y [ l\p) dp — > f(u) du 



J ./[0,oo)' i J/ ./[C^oo) 13 ° x 

for the first term of the left hand side of eq. (13.81) . □ 
Remark 3.9 The above two lemmata imply: 

lim Det [1 + v 7 ! - e-^Q^.Q^l - e^l = e K ^hf^ du . 

L— >oo 

Lemma 3.10 

lim Det [1 + y/l - e-f^Q K AL Q Vl - e~^(l + \fl - e-^Q x K kL Q^\ - e^)" 1 ] 

= Det [1 + \ffR\ff <8> COS 7TX 2 (cOS 7TX 2 , ® COS7r£ 3 (cOS7rX 3 , ■))] 

= Det[l+4- 1 v / 7^V / 7] • 

Here, the Fredholm determinant in the first member is for operators on L 2 (Al), the 
second determinant stays for operators on L 2 (I 3 ) and the third one for operators on 
L 2 (I) , see (TTTgp . 

Proof: Let U L : L 2 (A L ) -»■ L 2 (7 3 ) be unitary operator (U L g)(x) := L 2 g{L?x\, Lx 2 , Lx 3 ) 
for x G J 3 , g G L 2 (Al). Then integral kernels transform as 

{U L KUZ x ){x, y) = L^K{{L 2 x x , Lx 2 , La*), (LV, Lj/ 2 , Ly 3 )) , 

for x,y £ I 3 . It follows that 

L- 4 (U l Q K Al Q U l 1 )(x, y) = Q K A Q ((L 2 x u Lx 2 , Lx 3 ), (L 2 y u Ly 2 , Ly 3 )) 

8 sin(7rsi(2 _1 + xi) sin(7rsi(2 _1 + yi) cos(7rx 2 ) cos(7n/ 2 ) cos(7rx 3 ) cos(7n/ 3 ) 



r4 

si=i exp 



2m L 4 



(3.9) 

8m ^ cos[7rsi(xi - yi)} - cos[7rsi(l + x x + y x )\ . . . . . . . 

^ fl^2 2^ 2 , 2mL4A|^ op ~j COs(tTX 2 ) COS^) COs(7TX 3 ) COs(7T 

^ «i=l S l H 1 

= j/i) cos(7rx 2 ) cos(7ry 2 ) cos(7rx 3 ) cos(7n/ 3 ) 

uniformly in see (jl.74p and (lA.12j) . 
Then, we get the trace-norm limit: 

> y/fR\/J <E> C0S(7TX 2 )(C0S(7TX 2 ), (g) C0S(7TX 3 )(C0S(7TX 3 ) , •)) , 

which proves the lemma. □ 
Remark 13.91 and Lemma 13.101 yield Theorem ll.241 
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A Miscellaneous formulae 



In the appendix, the formulae used in the argument of this article are shown. Here, A 
is supposed to be a positive constant independent of L or Li,L 2 ,L 3 , while a positive 
constant B to depend on L or L\, L 2 , L 3 smoothly. We use the notation s — (si, s 2 , s 3 ) 
for sgN 3 and s — (si, s 2 ) for seN 2 , and so on. 

i° 



E 



1 



1 



L 2 

=(s 1 ,s 2 )eN 2 exp 
^(1,1) 



,4 



L 2 



B 



TV 

4A 



log(L 2 A B- 



0(1). 



(A.l) 



r,2 



L 2 

s=(s 1 ,s 2 )gn 2 exp 

st^(1 ,1) , ,S2'-odd 



A 



si-l , sS-l 



L 2 ^ L 2 



7T 

16A 



log(L 2 Afi- 1 ) + 0(1). (A.2) 



E ± 



L (L 2 B\ 

-rA—) +0{v - 



seN,s:octa,s>l exp 



.4 



s 2 -l 
L 2 



B 



(A.3) 



where 9? is defined in ( ll.70p . 
2° 

Suppose Li ^ L 2 ^ L 3 . The following convergence holds uniformly in x,y on any 
compact set as L\, L 2 , L 3 — > 00. 



E 



sew 3 exp 



.4 



-1 , si-l , s^-i 



+ 



B 



e 4A\p\ 2 +B(L=oo) _ l ■ 



(A.4) 



If Li ^ L 2 = L 3 . The following convergence holds uniformly in x,y on any compact set 
as Li, L 2 , L 3 — > 00. 



E 



s ,a(^)0 s ,a(2/) 



s6 m 3 exp 



5 



e 2Trip-(x-y)^p 
AA\p\ 2 +B(L=oo) _ I 



(A.5) 



Suppose L\ ^ L 2 ^ and 1 = (1, 1, 1). Then, we have 

E — 1 



s6N s L iM 3 exp 

S3* 1 



5 



- 1 



e .4|p| 2 +B(L=oo) _ l 



(A.6) 



as Li, L 2 , L 3 — > 00. 
3° 
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Suppose, Li ^ L 2 . Then 

T. - 



exp 



sf-l , si-i 



M + — 



^1 ^1 



5 



01 ^ A 1 



£2 L 9 V5 



OOog^AiT 1 / 2 )) 



(A.7) 



holds for large Li,L 2 , and L 3 . 

The following asymptotics also hold for large L: 



1 s\ + sl 



seN 2 exp 

1 /S? + S 2xl/4 



4 d-1 1 



0(l) + 0(L- 4 5 



L 2 V L 2 



^LL 2 



s=2 



4 



exp 



^-14 - 8 |-1 +g 



L 2 



0(1) + 0(L" 4 5 



exp 



A- 



s 2 -l 



L 2 



1) = 0(1) +0(L~ 5 / 2 B- 1 ] 





r s 2 - 1 




( ex p 




-0 



-1 



0(LAB" 1/2 ). 



(A.8) 
(A.9) 
(A.10) 

(A.11) 



Proof of flATTD : Put 



T 1 



L 2 

sgN 2 exp 



A 8 ;_i+ a 2_ 1 +g 



Because of 



L 2 



we get 

L.B. :-- 



1 \ 2 . /s, -l s,l 

for PiG (_, r 



dp 



(^N,j = l,2), 



1 



(o,oo )2 -(o,i/l] 2 exp[A(( Pl + 1/L) 2 + (ft, + 1/L) 2 - 2/L 2 ) + 5] - 1 

dp 



: U.B. . 



I? exp[3A/L 2 + B] - 1 y (0 , oo)2 _ (0 ,3/L] 2 exp[A(|p| 2 - 2/L 2 ) + 5] - 1 
For the lower bound (L.B.), we have 

* > / 

'(l/L,oo)2-(l/L,2/L]2 exp[A(H 2 - 2/L 2 ) + L?j - 1 J{\ p \ > 2V2/L, Pl ,p 2 >l/L} 



L.B. 



rn/2— sin 1 1/pL 

pc?j9 / d6 

2V2/L J sin- 1 1/pL 



1 



exp[A(p 2 - 2/L 2 ) + B] - 1 



12 



7T 

> - 



pdp 



psin l {l/pL)dp 



2 J2V2/L exp[A(p 2 - 2/L 2 ) + B] - 1 7 2v ^ /L exp[A(p 2 - 2/L 2 ) + B] - 1 " 
Here the first term is calculated as 

1 



7T 1 

o77T lo g 



2 2A ta 1 - e~ M / i2 - B 
For the second term, we use sin _1 x ^ tcx/2 to get 



£^4) + 0( i). 



the 2nd term > — 



7T 



dp 



7T 



LJ2V2/L exp[A(p 2 -2/L 2 ) + B]-l 
* 0(1). 



AL J 2V - 2/L p 2 -2/L 2 

For the upper bound (U.B.), we enlarge the integral domain and perform the angle 
integral to get 



U.B. < * 
2 



3/L 



exp[A(p 2 - 2/L 2 ) + B] - 1 L 2 exp[3A/L 2 + B\ - 1 



7T 



log 



AA 6 1 - e -7A/L 2 -s 
Proof of (JAJ3) : For 

1 



+ 0(1) 



7T 



log L A 



0(1). 



□ 



u- 



cxp 



^ (2n-l) 2 -l+(2m-l) 2 -l + g 



the inequalities 



2pi 



1\2 (2n-l) 5 



L 



L 2 



< 2p 



l\ 2 



2p 2 



1\2 (2m - I] 



L 2 



^ 2p 



L, 
1\ 2 



for pi G 



for p 2 G 



n — 1 n 



L L 

m — 1 m 



L ' L 



hold for n, m — 2, 3, • ■ • . The upper bounds of the above inequalities also hold for n = 1 
or m = 1. Thereby, we have 

dp 



(o,oo)Mo,W ex P[^(( 2 Pi + W + (2^2 + l/L) 2 - 2/L 2 ) + B] - 1 

00 00 

^ ^ 3^2,1 + 4 5*^1 + S^.m- 



(n,m)^(M) 



n=3 



n,m=3 



Calculation of these bounds can be performed in a similar way as in the proof of (lA.lj) . 
□ 
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Proof of (pO]) : Using 
it is obvious to get 



(A.12) 



El 



neN L exp 



L (I?B 
If 



A T \ A 



neN 



1 



cxp 



0(1). 



□ 



Proof of ( 1A.4|) : First, note that the inequality 

(2n - l) 2 - 1 (2m - l) 2 - 1 (2/) 2 - 1 



Li 



T 2 



T 2 

^3 



>Pl+p 2 2+Pt = \p\ 2 (A-13) 



holds for 



Pi e 



n — 1 n 



, P2 e 



m — 1 m 



L2 L2- 



P3 e 



Z- 1 z 



L3 L 3 



(n,m,l e N), 



and inequality 



(2n-l) 2 -l (2m-l) 2 -l (2/ — l) 2 — 1 2 2 2 , 2 

- + ji >p\+pI+pI= \P? (A.14) 



Li 



Lj 



holds for 



Pi e 



n — 1 n 



L\ 'La 



£>2 G 



m — 1 m 



L2 L2- 



P3 G 



/ - 1 Z 



L3 L 3 



{n,m,le M,Z^ 1), 



since La ^ L2 ^ L3. Now set 



. /(2n-l)vr (2n-l)7rsci\ . /(2n-lW (271-1)™ 
sin I 1 ; ) sm ' 



La 



'2rrm 2m7nr 2 \ (Ivrm 2mTry 2 
x sm I 1 sm 



Li 



•(2Z-1W (2Z-l)vrx 3 \ /(2/-1W (21 - 1W 
x sm I h r-^ I sm ' 



La 



La 



exp 



, (2rz-l) 2 -l (2m) 2 -l (2/- l) 2 - 1 



if 



Pi e 



n — 1 n 



La 'La 



, P2 e 



L| 



m — 1 mi 



L2 L2 



P3 e 



z - 1 z 



L3 L3 



(n,m,Z e N,Z ^ 1). 
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Otherwise, \I/£ eo (p; x, y) = 0. Similarly we can also define \^ e ^ e ^ e ^° etc. Put 

Hp) 



Z AW - 1 • 

Then, 

$ G L 1 ((0, oo) 3 ) and \^° L eo (p; x, y) | < $(p) 
hold. Due to the dominated convergence theorem, we get 

/ *HP;s,J/)dp— ► / r lim *r(P^,y)dp. (A.15) 

J(0,oo)3 </(0,oo) 3 

In the same way we also get the convergence of the integrals of \I/| eo (p; re, y), ty™ e (p; x, y) 
and so on, and finally 



s6 f« 3 G X P 



X4r + 4r + 4r)+ B 



- 1 

3 27rip-(x-j/) 



= 8 / {n ee (p; x, y)+*r(p; v)+- • -+*r(p; i/)) dp — ► / 4w^) _ i 

Here, note that the convergence in (1A.15I) is uniform in x, y on compact set. This follows 
from the following obvious observations: 

Suppose that F, Fi, F 2 , ■ ■ ■ G L 1 (IR 3 ) satisfy \\F n — F\\i — > 0, and that f(p, x), fi(j>, x), f(p, x), 
are bounded functions ofp,xE M. 3 satisfying \\f n — fW^ — > 0. 

Then, J F n (p)f n (p, x) dp —> f F(p)f(p, x) dp holds uniformly in x. □ 



Proof of (1A.5|) is almost the same as the proof of (1A.4|) 



Proof of ( 1A.6I) is the same as the proof of the limit of the integral for \P l in Lemma 
Proofs of the other formulae are straightforward. 
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